DAPNIA-04-37,ADP-04-03/T582 



Quark structure and nuclear effective forces 

P.A.M. Guichon^ and A.W. Thomas^ 

'SPhN-DAPNIA, CEA Saday, F91191 Gif sur Yvette, France 
'^Special Research Centre for the Subatomic Structure of Matter, University of Adelaide, SA 5005, Australia 

We formulate the quark meson coupling model as a many-body effective Hamiltonian. This leads 
naturally to the appearance of many-body forces. We investigate the zero range limit of the model 
and compare its Hartree-Fock Hamiltonian to that corresponding to the Skyrme effective force. By 
fixing the three parameters of the model to reproduce the binding and symmetry energy of nuclear 
matter, we find that it allows a very satisfactory interpretation of the Skyrme force. 
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The notion that quark degrees of freedom may play a 
role in low energy nuclear physics is largely unappreci- 
ated. The main reason is probably that the many body 
formulation of nuclear physics based on point like nucle- 
ons interacting through effective forces has proven quite 
successful. In this paper we take a radically different 
point of view, arguing that the nuclear effective force it- 
self is a direct manifestation of the quark structure of 
the nuclcon. To this end we formulate the quark meson 
coupling (QMC) model of the nucleus as a many body 
problem. This allows us to take the limit correspond- 
ing to a zero range force which can be compared to the 
Skyrme force 0. 

In the quark meson coupling model 0, 0] the essential 
step is to solve for the quark structure of the nucleon 
under the influence of the nuclear environment. For this 
one considers that, in a time averaged sense, a nucleus 
can be described as a collection of non overlapping quark 
bags representing the nucleons. (More recently the same 
ideas have been extended to a confined version of the 
NJL model The interactions of the quarks with the 
nuclear medium is represented by the exchange of mesons 
between the quarks of different nucleons, with coupling 
constants treated as free parameters. As explained in 
Ref.[3 the a field which is used here is not the chiral 
partner of the pion and the quark-u coupling does not 
break chiral invariance. 

As in our previous work 3] the scalar field is denoted 
(T(r), while is the time component of the vector field 
and both are taken to be time independent. In the nu- 
clear ground state, or for low energy excited states, the 
time dependence of the fields is driven by the Fermi mo- 
tion of the nucleons, so the typical frequencies are of the 
order of the Fermi energy, which can be neglected with re- 
spect to the high frequencies of the confined quark fields. 
Moreover, the space components of the u) field have their 
source in the velocity density of the nucleons, which is 
not a coherent quantity - in contrast to the nucleon den- 
sity p{r), the source of the time component. Therefore, 
when we solve for the nucleon structure under the influ- 
ence of the medium, the dominant effect comes from the 
time component 



Each bag is moving in the classical fields, a{r),uj{'r), to 
which the quarks are coupled. In the spirit of the Born- 
Oppenheimer approximation we solve the equations of 
motion of the quarks in a given bag for a fixed classical 
position, R(t), of its center. For this we use the known 
Lorentz character of the a and uj fields to transform to 
the instantaneous rest frame of the bag, where the static 
spherical cavity approximation is most appropriate. We 
then expand the fields around their values at the center 
of the bag, truncating the expansion at first order. The 
coupling of the quark to the constant part of the fields 
is solved exactly, because it amounts to a shift in the 
quark mass and energy. The remainder is treated as a 
perturbation. The rest frame energy-momentum is then 
transformed back to the nuclear rest frame with proper 
account of Thomas spin precession. Keeping only terms 
which are quadratic in the nucleon velocity, as we do 
systematically throughout this work, we find the follow- 
ing expression for the classical energy of a nucleon with 
position- momentum (i?, P) 
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with 5cjthe cj-nucleon coupling constant. The spin-orbit 
interaction, Vso, is defined below. To get the dynamical 
mass M*{R) one has to solve the bag equations in the 
field a{R). For our purpose it is sufficient to know that 
it is well approximated by the expression 

d / - x2 
2 

where (Af, g^) are the mass and the cr-nucleon coupling 
constant for the free nucleon and d = 0.22Rb - with 
Rb the bag radius. The last term, which represents the 
response of the nucleon to the applied scalar field (the 
"scalar polarizability" ) , is an essential element of the 
QMC model. From our numerical studies we know that 
the approximation |(2Jl is quite accurate up to g^a — 400 
MeV, which should be sufficient for our purposes. As a 
check we have also tried a cubic polynomial but found no 
significant changes in the results. 

The energy Q is the energy of one particular nucleon 
moving classically in the nuclear meson fields. Since, by 
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hypothesis, the bags do not overlap, the total energy of 
the system is the sum of the energy of each of the nucleons 
plus the energy carried by the fields. As the latter are 
static we can write: 
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with rricr, m^^ the masses of the mesons. 

To simplify the expression for Epf{R) we estimate the 
quantity ga-cr using the field equations 5Etot/ 5a{r) = 0. 
Neglecting the velocity dependent terms, setting M* « 
— QaO' and neglecting (Vcr)^ with respect to m^cr^, 
we find gaCT{r) ~ Gap'''{f), with Ga — gll'm'^ and the 
classical density is defined as p^'(r) — Yli S[f — Ri). 

From our previous studies in the Hartree approxima- 
tion, Ga ~ lOfm^, which yields gacr ~ 300 MeV at nu- 
clear matter density. Making a quadratic expansion of 
1/M* in powers of g^a and keeping the leading terms ^3 
we obtain: 
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If we define the scalar density as pf{r) = 



form 
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Ri), we can write the total energy in the 



drpf ( ga<y ~ ^{gacrf ) + I drp"^ g^u , (6) 



which will be our starting point for the many body for- 
mulation of the QMC model. We use the equations for 
the mesons, SEtot/ 5(j{r) = 5Etot/5uj{r) — 0, to ehminate 
the meson fields from the energy, leaving a system whose 
dynamics depends only on the nucleon coordinates. This 
is possible because the conjugate momenta of the static 
meson fields vanish identically. From Eq. ((Bj) we write the 
equations for the meson fields in the following forms: 



ga(T = Gapfil - dg^a) + \/'^g„cr/m 
g^Lo = Gujp"^ + V'^g^uj/m^ 
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where we have defined G^) — gt/'rn^- On the RHS of 
Eqs. (|7I8(I we have neglected the contribution of the func- 
tional derivative acting on the spin-orbit term in JHJ. 
This is because the latter was obtained as a first order 
perturbation and one can check that the resulting error 



in the final Hamiltonian is of higher order. If we insert 
the solutions Usoiir), <7soi{r) of Eqs. (|7I8|I in the expres- 
sion © for the energy we get, after some algebra, and 
omitting the irrelevant constant mass term: 

We do not attempt to use exact solutions of Eas. (|7l8|l 
as this would lead to an intricate many body problem 
that would be difficult to compare with standard nuclear 
physics approaches. Instead we first remark that, roughly 
speaking, the meson fields should follow the matter den- 
sity. Therefore the typical scale for the V operator act- 
ing on cr or cj is the thickness of the nuclear surface, that 
is about Ifm. In so far as lfm~^ <C (rn^,m^), which 
looks reasonable, we can consider the terms V^go-cr/m^ 
and V^^cjw/m^ as perturbations and, in these terms, re- 
place a and to by their first order approximation, that is: 
f/crCT w Gapf and g^i^ ~ G^jp'^K The next step in solv- 
ing for the a field is to solve Eq. Q iteratively , start- 
ing from the lowest order approximation, g^a = G^pf- 
When inserted in Eq. Q this series will generate N- 
body forces with convergence controlled by the parame- 
ter dg^a — 0.33 , according to our estimate. To simplify 
further we shall neglect the small difference between pf 
and p'^^ except in the leading term, the one which gen- 
erates the 2-body force. These approximations will not 
be difficult to improve but, as this leads inevitably to an 
effective interaction which is more complicated than the 
simple Skyrme force, we postpone this to future investi- 
gations. In summary the expressions we shall use for the 
field solutions are : 
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guOJsoi{r) = p" + GuP' 
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The rest of the derivation amounts to substituting 
gui^soi and gaCsoi into Eq. @ for the energy. As usual 
the density and the scalar density to some power contain 
infinite terms corresponding to the self-interaction of the 
nucleon. Formally these can be incorporated into a re- 
definition of the mass and coupling constants of the free 
nucleon. Since our model is devised to describe the modi- 
fication of the nucleon by the medium rather than the nu- 
cleon itself we simply remove them. This amounts to the 
replacements {J2iHr- Ri))'^ ^iy^j ^r-Ri)S{T^-Rj), 
which leads to the following many-body Hamiltonian, es- 
sentially equivalent to the QMC model: 
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Here Rij = Ri — Rj and Vi is the gradient with respect 
to Ri acting on the delta function. In Ea. H12|l we have 
dropped the contact interactions involving more than 4- 
bodies, because their matrix elements vanish for anti- 
symmetrized states. To shorten the equations we used 
the notation 5'^{ijk) for 6{Rij)6{Rjk) and analogously 
for 5^{ijkl) . For the spin orbit interaction we start from 
our previous result Q : 
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where W{R^) = M*{R,) + g^uj{R,){l - 2//,), = 0.9 
is the isoscalar magnetic moment and (Ji are the Pauli 
matrices. As this expression was derived as a first order 
approximation, it is consistent to evaluate it to the same 
order. So, on the RHS of Eq.lO we replace M* by 
M in the denominator, we approximate M* — M — g^a 
and we use the leading approximations for the meson 
fields - that is, gaCf = G^p^^, ff^w = G^p'^^ . As the spin 
orbit interaction is active only at the nuclear surface, the 
relevant density is rather small so these approximations 
are justified. 

The final step is to quantize the classical Hamilto- 
nian H12|) by making the replacement Pi —>■ — iV; . As 
usual we must deal with the ordering ambiguity which 
exists as soon as velocity dependent interactions are 
present. For the spin orbit interaction there is no am- 
biguity because all orderings give the same matrix ele- 
ments. The problem occurs only in the last term of the 
first line of Eq. H12II . There are 2 possible hermitian or- 
derings when Pi becomes an operator acting on the right: 
Ti = {PlS{R,j)+S{R^j)P^)/2 and T2 = P^S{R, - R,).P, 
. However, using integration by parts and the commuta- 
tion rules, one checks easily that the difference between 
the two orderings is of the form V-fS{Rij). Such an oper- 
ator is already present in Hqmc in the second and third 
lines of Eq. (|12|l and we see that choosing one ordering or 
the other is equivalent to a change of the meson masses. 
In practice we tried both orderings and checked that this 
is equivalent to a 50 MeV change of nia-, which is not 
so well known. Since, in any case, we intend to study 
the sensitivity of our results to m^, choosing Ti or T2 is 
immaterial and for definiteness and historical reasons we 
adopt the form T2. 

To complete our effective Hamiltonian we now include 
the effect of the isovector p meson, which can be done 



by analogy with the w meson. If we let b" {a = 1,2,3) 
be the time component of the field and r" the isospin 
Pauli matrices, then the only changes are the replacement 
g^ui{R) — > g^uj{R) + gph{R).f/2 in the expression ^ for 
the nucleon energy and gujijj(R)(l — 2ps) 9i^^{R)i^ — 
2ps) +5p(l — 2pv)b{R).f/2 in the expression (fO|l for the 
spin orbit interaction, with gp the free p-nucleon coupling 
constant and py — 4.7 the nucleon isovector magnetic 
moment. If we define Gp — gp/rn^p with rrip the mass of 
the p meson, our quantum effective Hamiltonian finally 
takes the form : 
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with Aij =Gcr + {2ps - l)Gu + {2pv - l)GpT,.fj/4. 

To fix the free parameters of the model, that is G^, G^ 
and Gp, we have computed, using the Hamiltonian H14|l . 
the volume and symmetry coefficients of the binding en- 
ergy per nucleon of infinite nuclear matter : Eb/A = 
ai + ai[N ~ ZY jA? . We have used the experimental val- 
ues a\ = — 15.85MeV, = 30MeV and the saturation 
condition daxjdp^po) = with po — 0.16fm~'^. In order 
to avoid the proliferation of tables we show only the re- 
sults corresponding to the bag radius Rb = 0.8 fm, which 
is realistic. We have used the physical masses, = 
782MeV, m,p = 770MeV and we allow ma to take the val- 
ues SOOMeV and 600MeV, which is a commonly accepted 
range. We get, in fm^ G^ = 12.63, G„ = 9.62, Gp = 9.68 



for 



= 500MeV and G„ = 11.97, G<^ = 8.1, G„ = 6.46 



for vria ~ 600MeV. These values are larger than in the 
Hartree approximation because the exchange terms tend 
to cancel the direct terms of the matrix elements, thereby 
forcing larger couplings to fit the data. 

As a practical test of the capacity of our model to in- 
terpret a large body of nuclear data, we compare it with 
the effective Skyrme interaction. Since, in our formula- 
tion, the medium effects are summarized in the 3 and 4 
body forces, we consider Skyrme forces of the same type, 
that is without density dependent interactions. They are 
defined by a potential energy of the form : 

Y^t^Y. SiR,j)6iRjk) + toY.il + xoPa)S{R,j) 
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with Vij — Vi — Vj. There is no 4-body force in 
Eq. H15|l and we shaU show its possible impact by set- 
ting its strength equal to zero in Hqmc- Since the spin 
exchange operator, P^, mV multiplies a contact interac- 
tion, its action on an antisymmetric state it is equivalent 
to minus the isospin exchange operator. Comparison of 
Eq. (^ni) with the QMC Hamiltonian (^H) allows one to 
identify 

to = -G, + G^-^, h^MGl, (16) 

For the other parameters we cannot make a direct iden- 
tification because, as our effective Hamiltonian is derived 
in the rest frame of the nucleus, its momentum depen- 
dent pieces violate Galilean invariance. This is irrele- 
vant since it is devised for variational calculations where 
the trial state also violates Galilean invariance but to 
make the identification we need to compare the respective 
Hartree-Fock Hamiltonians rather than the interactions 
themselves. To this end we make some simplifying as- 
sumptions which do not significantly damage the physics 
but avoid unnecessary technical complications. First, we 
restrict our considerations to doubly closed shell nuclei 
with N = Z. Second we assume that one can neglect the 
difference between the radial wave functions of the single 
particle states with j = 1 + 1/2 and j = I — 1/2 . This 
amounts to treating the spin orbit interaction to first or- 
der, which is sufficient for our purposes. By comparing 
the Hartree-Fock Hamiltonian obtained from Hq mc a-nd 
that of Ref. corresponding to the Skyrme force , we 
obtain the relations : 

3^i + 5^2 = S+4(%-%)+3%, (17) 
- 9^1 = + 28 f % - - 3% , (18) 

Wo = (^5G, + 5(2/i. - 1)G. + J(2m. - l)Gp) . 

(19) 



TABLE I: QMC predictions compared with the Skyrme force. 





QMC 


QMC 


Skill 


QMC(N=3) 


m<,(MeV) 


500 


600 




600 


to(MeVfm^) 


-1071 


-1082 


-1129 


-1047 


Xo 


0.89 


0.59 


0.45 


0.61 




16620 


14926 


14000 


12513 


3fi -h5f2(MeVfm^) 


192 


475 


710 


451 


5t2 -9fi(MeVfm^) 


-7622 


-4330 


-4030 


-4036 


Wo{MeVim^') 


118 


97 


120 


91 


K{MeV) 


327 


327 


355 


364 



In Tabled we compare our results with the parameters 
of the force Skill , which is considered a good represen- 
tative of density independent effective interactions. We 
postpone a more extensive comparison to future work. 
Instead of ti,t2, we show the combinations 3ti + 5^2, 
which controls the effective mass, and 5^ — 9ii, which 
controls the shape of the nuclear surface Even bear- 
ing in mind that, since we use the parameters (ai, 04, po) 
as input, not all the numbers in Table are predic- 
tions, we see that the level of agreement with Skill, for 
Too- = 600MeV, is still impressive. An important point is 
that the spin orbit strength, Wq, comes out with approxi- 
mately the correct value, independent of the a mass. The 
last column (N=3) shows our results when we switch off 
the 4-body force. The main change is a decrease of the 
predicted 3 body force. Clearly this mocks up the effect 
of the attractive 4-body force which may then appear less 
important. However, this is misleading if we look at the 
compressibility of nuclear matter, K — dp^d^ai/dp^ (last 
line of Table P), which decreases by as much as 37MeV 
when we restore this 4-body force. The value we find, 
K = 327MeV, is still a little too large with respect to the 
experimental range (200 + 300MeV) , but several simpli- 
fications made for this presentation can be eliminated in 
future work. Moreover we have not yet included the long- 
range force of the pion, which will reduce K by about 
20MeV according to a preliminary calculation. This too 
will be investigated in future workj^. 

In summary, we have demonstrated a remarkable 
agreement between the phenomenologically successful 
Skyrme force. Skill, and the effective interaction cor- 
responding to the quark meson coupling model - a re- 
sult which suggests that the response of nucleon internal 
structure to the nuclear medium does indeed play a vital 
in nuclear structure. 
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